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Some fundamental and formal aspects of the quantum dwell time are reviewed, examples for free 
motion and scattering off a potential barrier are provided, as well as extensions of the concept. We 
also examine the connection between the dwell time of a quantum particle in a region of space and 
flux-fiux correlations at the boundaries, as well as operational approaches and approximations to 
measure the flux-flux correlation function and thus the second moment of the dwell time, which is 
shown to be characteristically quantum, and larger than the corresponding classical moment even 
for freely moving particles. 
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Time observables in quantum mechanics have a long and debated history (35j. In spite of the fact that random 
time variables, measured after a system is prepared, are common in laboratories, most often it has been argued that 
questions about time in quantum mechanics should best be left alone, as illustrated by the frequent reference to Pauli's 
i— i. theorem. Alternatively, the emphasis has been laid on characteristic times, i.e. single time quantities characterizing a 
' process such as tunneling, or decay. This, in many ways, runs counter to the usual procedure in quantum mechanics, 
where additionally to the average value of a quantity we require prediction of higher order moments of that quantity; 
in other words, the probability distribution. 

With regard to the timc-of-arrival observable several such distributions have been proposed and studied, see volume 
1 of "Time in Quantum Mechanics" fl2l |. or [3C| . In this chapter we analyze the dwell time, which appears to be a 
Q-f much simpler time observable because the associated operator will indeed be self-adjoint (over the adequate domain). 
— — At first sight, it could be thought that this statement contradicts Pauli's theorem, which asserts that no self-adjoint 
time observable can exist with canonical commutation relations with a semi bound Hamiltonian. However, the dwell 
time is an interval quantity, as opposed to the instant quantity that the time of arrival describes. Its associated 
I 1 operator, therefore, should commute with the Hamiltonian, as opposed to presenting canonical commutation with it. 
| The dwell time of a particle in a region of space and its close relative, the delay time (45j . are rather fundamental 
quantities that characterize the duration of collision processes, the lifetime of unstable systems [13j . the response 
to perturbations [2(j, ac-conductance in mesoscopic conductors 0, or the properties of chaotic scattering [3l|. In 
y—t • addition, the importance of dwell and delay times is underlined by their relation to the density of states, and to the 
virial expansion in statistical mechanics [8]. We could thus hardly fail to study and characterize in detail such a 
ON , prominent quantity. For a sample of theoretical studies on the quantum dwell time see [!, [H, H, [ll|, US HE H3, HE 
HI, EE EE SI, EE US HH H3] ■ A recurrent topic has been its role and decomposition in tunneling collisions. Instead, 
we shall focus here on a different, so far overlooked but fundamental aspect, namely, the measurability and physical 
implications of its distribution, and its second moment. 

Despite the nice properties of the dwell-time operator, the relevance of the concept and average value in many 
different fields, or the apparent formal simplicity stated above, the dwell time is actually rather subtle and remains 
elusive and challenging in many ways. In particular, a direct and sufficiently non-invasive measurement, so that 
the statistical moments are produced by averaging over individual dwell-time values, is yet to be discovered. If the 
particle is detected (and thus localized) at the entrance of the region of interest, its wavefunction is severely modified 
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("collapsed"), so that the times elapsed until a further detection when it leaves the region do not reproduce the ideal 
dwell-time operator distribution, and depend on the details of the localization method. Proposals for operational, i.e., 
measurement-based approaches to traversal times based on model detectors which study the effect of localization have 
been discussed by Palao et al. [3!| and by Ruschhaupt [4l[ . For attempts to measure the dwell time with continuous or 
kicked "clocks" coupled to the particle's presence in the region see (ll.l49j. All operational approaches to the quantum 
dwell time known so far have provided only its average, and indirectly by deducing it from its theoretical relation 
to some other observable with measurable average. The average is obtained for example by a "Larmor clock", using 
a weak homogeneous magnetic field in the region D and the amount of spin rotations of an incident spin-i particle 
An optical analogue is provided by the "Rabi clock" [4<] . It can also be deduced from average passage times 
at the region boundaries [29| , as well as by measuring the total absorption if a weak complex absorbing potential acts 
in the region [ljl [H, HH . This last setup could be implemented with cold atoms and lasers as described in [ll|, EH 
and will be discussed in Sect. IIVB1 

The chapter is organized as follows: the first sections are devoted to fundamental and formal aspects (Sect. 2), 
examples (Sects. 3 and 4), or extensions (Sect. 5) of the dwell-time concept and operator, whereas Sect. 6 tackles the 
relationship between the moments of the dwell-time operator and flux- flux correlation functions (ffcf) [371 ] . generalizing 
an approach by Pollak and Miller |4(J. They showed that the average stationary dwell time agrees with the first 
moment of a microcanonical ffcf. We shall see that this relation holds also for the second moment, but not for higher 
moments, and extend their analysis to the time-dependent (wavepacket) case. We shall also discuss a possible scheme 
to measure ffcf 's, thus paving the way towards experimental access to quantum features of the dwell-time distribution. 

II. THE DWELL TIME OPERATOR 

Unlike other time quantities, there has been a broad consensus on the operator representation of the dwell time 
[l3L I20I ]. For one particle evolving unitarily with Hamiltonian H in region D, which we limit here to one dimension 
for simplicity, D = {x : x\ < x < X2}, it takes the form 

/oo />oo 
dtxn(t)= dte im / h XD (x)e- iHt / h , (1) 
-00 J — 00 

where xr>W is the (Heisenberg) projector onto D, and xd{%) = Xd(0) = dx\x)(x\. Without delving too far in the 

functional analysis definition, i.e. into the proper description of its domain and of its adjoint, we can see that To will 
be self-adjoint (we shall come back to this issue after studying the specific example of the free particle Hamiltonian 
in Sect. IIII[) . At any rate, it is clear that, at least formally, this operator commutes with the Hamiltonian, as can be 
seen from what follows: 1 

/oo , 
dTC iHT / h X D(x)e- iHiT+t V h (2) 
-OO 

/OO 
dTe m-r-t)/n XD ^ )e -iH T /n = e -im/nf D 
-00 

The commutation of To and the Hamiltonian leads us to search for the eigenfunctions of dwell time in the stationary 
eigenspaces of the latter. Let a be the degeneracy index for these eigenspaces, such that H\E,a) = E\E,a). We 
easily obtain the matrix elements of To in the corresponding eigenspace, 

T D \E,a)=2TchJ2(E,P\XD@)\E,a)\E,(3) . (3) 


We may thus reduce the problem of eigenvalues and eigenfunctions of the dwell-time operator to a set of matrix 
diagonalization problems in each of the eigenspaces of the Hamiltonian. 

Except in Sect. |V]we shall assume that the Hamiltonian holds a purely continuous spectrum with degenerate (delta- 
normalized) scattering eigenfunctions \4>±k) corresponding to incident plane waves | ± k), with energy E = k 2 Ti 2 / (2m), 
normalized as (k\k') — (4>k\(l>k'} = 5(k — k'). 



Time-limited versions of the dwell-time operator such as J °° dt e lHt / h XD (5) e ' Ht / n do not generally commute with H, see [13| | or Sect. 

ivbi 



Following the same manipulation done for the S'-operator in one-dimensional scattering theory [29j | , it is convenient 
to define an on-the-energy-shell 2x2 dwell-time matrix T, by factoring out an energy delta, 



m 



where E' = k /2 hy(2m) and 



Tfc/t' = {<f>k\XD(x)\<f>k> 
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"im' 



E = E' 



(4) 



(5) 



In particular, T k k is the average dwell time for a finite space region defined by Biittiker in the stationary regime 

1 



Tfefe 



\m\ 



dx \4>k{x)Y 



(6) 



where j(k) is the incoming flux associated with \4>k)- 

An intriguing peculiarity of the quantum dwell time is that the diagonalization of T at a given energy provides in 
general two distinct eigenvalues t±(k), k > 0, and corresponding eigenvectors \t±(k)}, even in cases in which only a 
single classical time exists, such as free motion, or transmission above the barrier; some explicit examples are discussed 
below. A consequence is a broader variance of the quantum dwell-time distribution compared to the classical one. 

The quantum dwell-time distribution for a state \ip) = \ip(t = 0)), is formally given by 



n^(r) = (tP\6(T d - r)|V) , 



(7) 



since the self-adjointness of the dwell-time operator implies that the spectral theorem applies, and that operator 
moments coincide with the moments of the distribution. We shall not consider in this work other sources of fluctuations 
such as mixed states or ensembles of Hamiltonians. In these two cases one could compute distributions of average 
dwell times, whereas here we shall be interested in the distribution of the dwell time itself, but only for pure states 
and a single Hamiltonian. 

On computing the distribution of dwell times we run into the difficulty, mentioned above, that the dwell-time 
operator is multiply degenerate, and that it is not altogether easy to de-scramble the relation between a given t 
and the corresponding set of momenta. It is useful in this case to profit from the spectral theorem to compute the 
generating function of the dwell-time distribution, defined as 



Normalizing the dwell-time eigenvectors so as to have the resolution of the identity 
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can be written as 
and, hence, 



dk 



^\ t+{k)) ( t+ (kM + e tot - « (#_ (A:)) (t_ (k)\^) 



(8) 



(9) 



(10) 
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dk [5(t~t + )Mt + )\ 2 + 5(t-t_)Mt_)\ 2 ] 



(11) 

The average of the dwell time for a wavepacket can be given in terms of the position probability density in cor- 
respondence (unlike its second moment, as shown below) to the expression for an ensemble of classical particles. It 
reads M,M 



dt 



dx\il)(x,t){' 



dfc|(£# m )| 2 T fefc , 



(12) 




FIG. 1: Eigenvalues of the dwell-time-operator as a function of q = kL, in units mL 2 /h, for a freely moving particle in a region 
of length L 

where ip(x,t) = dk (4>k\ip) ex P( — ifik 2 t/2m)(j)k{x) is the time-dependent wave packet and we assume, here and in 
the rest of the Chapter, incident wavepackets with positive momentum components. To write Eq. (fTZ)) use has been 
made of the standard scattering relation {4>k\ip) — (k\ip m ), where (x\k) = {2ir)~ 1 / 2 exp(ifoc), and ip m is the freely- 
moving asymptotic incoming state of tp. Space-time integrals of the form (|12p had been used to define time delays by 
comparing the free motion to that with a scattering center and taking the limit of infinite volume [l5j |. 

The second moment takes the form, as before for wavepackets with incident positive momentum components, 



n) = 



dk(T 2 ) kk \{kW 



dk- 



47r 2 m 2 
h 2 k 2 



dfc(|T fcfc | 2 + |T fc _ fc | 2 )|(£# m }| : 



( r dx\^)\- 



dx (f)l(x)<f>-k(x) 



\{kW n )V, 



(13) 



with a term without classical counterpart. 



III. THE FREE PARTICLE CASE 

In order to get a better grasp of the abstract results, it is adequate to illustrate them with the explicitly solvable 
free-particle case for a region D that extends from x\ — to xi = L. Indeed the humble freely moving particle turns 
out to be rather interesting and surprisingly complex with regard to the dwell time. The free-particle Hamiltonian 
is doubly degenerate; we can choose the degeneracy indices to coincide with the sign of the momentum, and, on 
computation, with due attention to the different normalization of the energy and the momentum eigenfunctions, we 
find the following generalized eigenfunctions of the free dwell-time operator: 



\t±{k),D) 



1 



V2 



[|*)± 



AkLi 



k)] 



(14) 



where 



mL ( 1 



(15) 



are the corresponding eigenvalues, in clear contrast to the classical time t c \ ass = mL / (\k\Ti) , see Fig [TJ It is important 
to notice that the inverse function is multivalued. Due to this multivaluedness we have kept generalized eigenfunc- 
tions with the dimensionality of |fc), instead of normalizing them to the delta function in dwell times. With this 
normalization it is easy to check that the resolution of the identity has the form Q . 

Interestingly, i_(fc) tends to zero as k — > 0, one more very non-classical effect that is better understood from the 
coordinate representation of the eigenvectors, 



(x\t+) 
(x|t_) 



e ife_L/2 

TT 1 / 2 
i e ifei/2 

TT 1 / 2 



cos[k(x - L/2)] , 
im[k{x - L/2)] 



(16) 



symmetric and antisymmetric with respect to the interval center. (a:|t_) vanishes at L/2, so the particle density tends 
to vanish in D with longer wavelengths as k — > 0. 

The dwell-time operator for an interval D' — [a, a + L] is obtained from the one above by translation, 



T D . =e- iak T D e iak , (17) 
and, as a consequence, the eigenvalues are not modified, while the new eigenfunctions are easily computed as 



\t±{k),D') = e- iak \t ± (k),D 
1 

" V2 



— iak 



\k) ±e ifc ( a+L >| -k) 



It is straightforward to compute directly the action of in the wavenumber representation, 



\TdW) = 



mL 
\k\h 



f/>(k) + e' lkL —sm (kL) ip(-k) 
kL 



(18) 



(19) 



where ip(k) = (k\ip). This allows us to study the functional aspects of the operator. In particular, it is easy to check, 
using the requirement of normalizability, that the domain of the operator is given by functions such that ip(k)/k — > 
as k — > 0. The requirement of symmetry does not add further limitations to the domain. As to the deficiency indices, 
they are computed to be (0, 0). These computations are carried out for the interval D = [0, L], but, given the unitary 
equivalence of other intervals of the same length, see Eq. (fTT|) . the same results carry over for regions composed of 
an arbitrary number of closed intervals. 

It should come as no surprise that functions in the domain of the (free particle) dwell-time operator must vanish at 
p = fast enough, since the characteristic evolution of a generic wavefunction is dictated by the free propagator, with 
its £ -1 / 2 temporal behavior. This entails the divergence of the dwell time unless the state has no p = component 
and the amplitude decay is faster (for a simple analysis, see [13]). The divergence of the average dwell time for states 
with nonvanishing p — components also occurs for ensembles of classical particles. 

To calculate the distribution II,/, (i), see Eqs. (TTU|) and (fTTj) , we have in this case the characteristic function 



dke i u mL/(\k\h) cos 

X) 

dke iw n L/(\k\h) gin 



ijjm 
Lum 



sin(fcL) \i/>(k)\ 



\k\kh 



s'm(kL) 



e~ ikL ij(k)4>(-k) 



(20) 



The distribution EL^, (t) has support only on the positive semi-axis, as can be seen from the explicit expression for the 
eigenvalues. This, in turn, is a nontrivial check of the correctness of the definition. 

A further initially surprising result of this analysis is that for highly monochromatic wavepackets (i.e., highly 
concentrated in one point in the momentum representation), the probability density for dwell times is generically 
bimodal because of the two eigenvalues t±(k) expounded in Eq. (fT^j) . An experimental verification of the eigenvalues, 
and of the quantum nature of the dwell time, could be realized with the aid of a matter-wave mirror located at a 
point X > L, reflecting an incident plane wave \k). The resulting standing wave would take the form, up to a global 
phase factor, 



|V.x) = |fc)-c— |-fc) 
We may now compute the average dwell time between and L, 

(1>x\T\il>x) 



(21) 



T kk +T_ fc _ fc - 2Re(e- 2lkX T k _ k ) 

J 1 - en 

kin 



kL 



(22) 



which oscillates between the maximum and minimum values 2t±, see (|15[) : they occur at specific locations of the 
mirror, namely, 



(n integer such that X > L), for which the standing wave \ipx) becomes proportional to \t±). The proportionality 
factor 2 1 / 2 accounts for the fact that the extrema correspond to twice the eigenvalues, which is easy to interpret 
physically with reference to a corresponding classical scenario: the classical particle under similar circumstances 
would traverse the region twice, first rightwards and then leftwards. For any X between the priviledged values X± 
given above, \ipx) is a linear superposition of the dwell time eigenvectors and thus the resulting average dwell time 
lies in a continuum between two times the eigenvalues. In a proposed experiment a highly monochromatic continuous 
beam would be sent towards the mirror and the particle density could be measured between and L by fluorescence 
or other means. The oscillations of the signal as a function of X would be in sharp contrast to the classical case, for 
which the the beam density and dwell time would remain unaffected by a change of the mirror's position. 



A. A comparison with classicality 



We have already noticed some of the peculiarities of the quantum dwell time compared to the classical dwell time. 
Here we elaborate this comparison further. Consider a classical ensemble of free particles, described by the initial 
probability density on phase space, F(x,p). The dwell-time distribution for this case is given by 



n c iass(i) = J dxdpS f* _ ~^j~) F ( x > 



P) 



mL 



dx 



mL\ „ / mL 



(25) 



that is to say, the marginal momentum distribution evaluated at mL/t and —mL/t, and multiplied by the normaliza- 
tion factor mL/t 2 . On the other hand, the distribution (jTTJ) obtained from Eq. ([2TH includes effects of interferences 
between positive and negative momentum components, in two different ways. In the first place, there is the obvious 
interference of the last line of (f2"D|) ; but, additional to this, the argument of the cosine also reveals these effects. In 
order to see better this point, let us examine a different operator, 



t D := X + T D X + + A_T D A_ = mL/\p\ 



(26) 



where X± are the projectors onto the positive/negative momentum subspaces. Notice that X± do not commute with 
Td. The last form in Eq. (f2"6"|) . specific of free motion, is particularly transparent and reproduces the one of the 
classical dwell time. The eigenfunctions are | ± k), k > 0, and the corresponding eigenvalues are twofold degenerate 
and equal to the classical time, mL/Ti\k\. The distribution of dwell times for this operator, and for positive-momentum 
states, is given by 



mL 



(27) 



which coincides with the classical distribution for initial wave functions whose support is limited to positive momenta. 
In this manner we see that the distribution incorporates interferences even if only positive momenta are present. 
To be even more explicit, observe that, for a state ip with only positive momenta components, the second moment of 
the dwell-time operator {1} is given by 



dfc(T 2 ) fefc |^(fc)| 5 



dfc(|T fcfc | 2 + |T fc _ fc | 2 )|-0(fc)| ; 



2 t1 



dk 



m L 

k 2 n 2 



k 2 L 2 



sin 2 (kL) 



(28) 



which contrasts with 



o 



00 „ m 2 L 2 
dk—^- 



k 2 h l 



thus indicating that indeed the second term is due to quantum interference. 

As indicated above, there is another rather striking way of seeing this point, by considering a highly monochromatic 
wavepacket. The classical distribution would have a very sharply defined peak around mL/po, where po is the central 
momentum of the wavefunction, po = kgh; on the other hand, the quantum distribution would have two sharp peaks, 



centered on t+(ka) and i_(fco) respectively. The distance between peaks goes to zero as pa increases, as was to be 
expected since that is the classical limit. 

The on-the-energy-shell version of to, t, is also worth examining. By factoring out an energy delta function as in 
Eq. (ID) we get for a plane wave \k) the average tkk = mL/(hk), which is equal to Tkk, but the second moment differs, 
(t 2 )fefc = (tkk) 2 = (Tkk) 2 < (T 2 )/c/c, see Fig. 1; in other words, the variance on the energy shell is zero since only one 
eigenvalue is possible for t. Contrast this with the extra term in Eq. (|28p . which again emphasizes the non-classicality 
of the dwell-time operator To and its quantum fluctuation. 



IV. THE SCATTERING CASE 



Let us now study the dwell time in a potential barrier or well without bound states, that is, the dwell time in an 
interval which coincides with the finite support of the potential V(x), which presents no bound states. For simplicity, 
we shall assume as before that this interval starts at point x = and has length L. We shall use the complete basis 
of incoming scattering stationary states, \k + ), with the following position representation, for k > 0, 



whereas, for k < 0, we have 



(x\k+) = 



(x\k + ) 



1 



Akx 



R l (k)e 



-ikx 



ikx 



2^ 1 T l (k)e 



1 / T r (-k)e [kx 



for x < 
for x > L ' 



'2-k 



Akx 



R r (-k)e 



-ikx 



for x < 
for x > L 



(29) 



(30) 



(the superscripts I, r in T , R l , and T r , R r stand for left and right incidence). We have omitted any explicit expression 
in the interval [0, L] because of the potential dependence. The eigenvalues and eigenstates of the dwell-time operator 
can be formally computed using H \k+) = (k 2 h 2 /2m) \k + ). To that end define 



£(*) = 



(k+\ X o(x)\k+) 



where 



and 



2a(k) 

a(k) = \(-k+\ X D(x)\k+)\ 
-k+\ XD (x)\k+) 



j v (k) 



a(k) 



Additionally, for the sake of compactness in later formulae, let 

^) = \[i{-k)-m\ ■ 

We then have that the eigenvalues of the dwell-time operator are given by 



t±(k) 



2-Kma(k) 
\k\Ti 



±y/l + ^(k) 



while the eigenfunctions are 

\t ± (k)) = N ± [\k+) + e iv(k *> [fi{k)± ^1 + M 2 (fc)] 
and the normalization is given by 



1 + /i 2 ± /iv/l + /i 2 (fc) 



-1/2 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



In fact, we can relate the quantities cr(fc) and tp(k) to scattering data whenever the support of the scattering 

potential is completely included in the region D. Let Xd(x) be the complementary function to xd(x). We can compute 



(k' + \x(x) \k + ) explicitly, and, using unitarity and the conditions this imposes on the scattering amplitudes, we are 
led to the following expressions (where we omit the arguments of the scattering amplitudes, since all are evaluated at 
k, and denote derivative with respect to k as d k ): 



(k+\ X D(x)\k+) 



L 



2n 
L r 



— [R l d k R l + T l d k f l ] + -r—r [R l — R l ] , 
2n 1 J Airk 1 J 



(-k+ \ X d(x)\ -k+) = — [l + \R r \ 2 ] + -L [R r d k R r + T r 8 k f r ] 



Airk 



-2ikL R r _ e 2ifci^rj 

i 



(-k+ \ XD (x)\ k+) = ±-T l R r + -L [R l 8 k T r + T l d k R r ] 



ink 



[T r 



Notice that the second term of each of the previous expressions can be identified with the corresponding on-shell 
matrix elements of —(\m/kT%)S^d k S. That is, with Smith's delay time! 



A. The Square Barrier 

For the specific case of a square barrier, in which V(x) — Vq[8(x — L) — 9(x)], with 6(x) being Heaviside's unit step 
function, we can put to good use the symmetries of the Hamiltonian, namely x — > L — x and k — > —k. These are 
realized on the scattering states as 



(L — x\k + ) 



AkL 



whence £(fc) = £(— k), and = 0. Furthermore, 



(x\-k+) 



(k+\ XD (x)\k+) 



L\T{k)Y 
2ttk 2 



k 2 - 



mVv, 



1 H sin (2kL) 

2kL 



where k = y k 2 — 2mVo/h 2 , real and positive for positive above-the-barrier momenta, and with positive imaginary 
part for tunneling momenta. After some cumbersome algebra, one readily obtains 



t±(k) = 



\T{k)\ 2 mL 
2fi\k\K 2 



[k 2 + k 2 ± (k 2 - k 2 ) cos kL] 



1 ± — sin kL 
kL 



which is to be compared with (j 1 5 j) ; namely, the free particle result gets modulated by the transmission amplitude and 
a potential dependent oscillatory factor. For high energies the modulating term tends to one, and we recover the free 
particle case, as one should. For small momenta, both eigenvalues tend linearly to zero with the momentum, unlike 
the free case where one of them diverges. It is also relevant that the dwell-time eigenvalues are bounded above, and 
therefore the average value is also bounded! In order to have a better grasp of this result, it is convenient to express 
the eigenvalues in dimensionless terms (k = q/L, Vq = h 2 Q 2 /(2mL 2 ) ): 



t±{k) 



2mL 2 



q± 



\/q 2 -Q 



= sin \/q 2 - Q 2 



ft 2q 2 -Q 2 ±Q 2 cos vV - Q 2 ' 



see Fig. H 



B. Average dwell time from fluorescence measurements 

In this subsection we shall model the measurement of the average dwell time t& = (To) of an ultracold atom at a 
square barrier created by a laser shining perpendicularly to its motion with homogeneous intensity between and L. 

Consider a two level system coupled in a spatial region to an off-resonance laser with large detuning, A 3> 7, O, 
where A is defined as the laser frequency minus the frequency of the atomic transition, 7 is the decay constant (inverse 
lifetime or Einstein's coefficient), and is the Rabi frequency. 




q 



FIG. 2: Dwell time eigenvalues for the square barrier, in units of mL 2 /%, for Q = 8 



The amplitude for the atomic ground state up to the first photon detection is governed then by the following 
effective potential, see [27],[38l or Chapter 

w <w 

so that the average detection delay (life time of the ground state if the atom at rest is put in the laser- illuminated region) 
is, see e.g. [36| . 4A 2 /fl 2 j. Whereas 7 is fixed for the atomic transition, Q and A may be controlled experimentally, 
and the ratio f2 2 /A can always be chosen so that the real part of V remains constant. This still leaves some freedom 
to fix their exact values which we may use to set the imaginary part. If we do so making sure that at most one 
fluorescence photon is emitted per atom, i.e., td ^> 4A 2 /f2 2 7, so that the fluorescence signal produced by an atomic 
ensemble will be proportional to the absorption probability A, this signal provides, after calibration to take into the 
detector solid angle and efficiency, and approximation for the derivative (|40|) and therefore to the average dwell time 
for the potential (|3"5|) as 

t d » HA/(2Vj) . (39) 

This follows by integrating —dN/dt = (2Vi /h)(ip(t)\xD(x)\ip(t)) over time, N being the surviving norm, and A = l — N 
the absorption (fraction of atoms detected). In the limit Vi — > Q and for highly monochromatic incidence, the average 
(stationary) dwell time at the real potential is obtained, 

T kk = lim (h/2)d Vl A(k) , (40) 
vi — >o 

where A(k) is the total absorption probability for incident wavenumber k. The equivalence of this quantity with 
(t+(k) + t-(k))/2 can readily be checked. 

One may think of relaxing the one-photon condition to get a proxy for the dwell time of an individual atom of the 
ensemble from the photons detected in an idealized one-atom-at-a-time experiment. For Vr negligible versus E, it 
could be expected that for some regime this distribution of emitted photons would also be bimodal. For the bimodality 
to be observed, the characteristic interval between modes (h/E, where E is the particle's energy) should be greater 
than the characteristic interval between successive emission of fluorescence photons, but these conditions and A > 7 
are not compatible, and similar difficulties are found for on-resonance excitation. A pending task is the application 
of (deconvolution or operator normalization) techniques which have been successfully applied to the arrival time, at 
least in theory. 

Figure [3] shows the exact dwell time and approximations for several values of Vj calculated for a transition of Cs 
atoms (the details are in the figure caption). A larger Vi implies larger errors but also a stronger signal. In practice 
the minimal signal requirements will determine the accuracy with which the dwell time can be measured. Figure 0] 
shows the relative error of the dwell-time maxima versus the corresponding absorption probability. In this figures the 
beam is monochromatic. We can check the reality of the quantum prediction at hand, differing from the classical one; 
namely, that for all ingoing waves the quantum mechanical dwell time is bounded, unlike the classical one. 
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FIG. 3: Exact average dwell time (solid line) for Cs atoms crossing a square barrier of width 2/^m and height 8.2674 x 10 3 s -1 7i 
(in velocity units, 0,28 cm/s) versus incident velocity. Approximations are calculated with Eq. (I39|l for Vi = 3.307s _1 fi = Vi 
(indistinguishable from the exact result, A = 25007, fl = I.577), 10 Vi (double dotted-dashed line, A = 2507, Q = O.57), 
10 2 Vi (dashed line, A = 257, Q = .I67), and 10 3 Vi (dotted dashed line, A = 2.57, ^ = O.O57). The transition is at 852 rim 
with 7=33.3 xl0 6 s _1 ; A and Q are obtained from Eq. fl38j| 
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FIG. 4: Relative errors, calculated from the maxima of exact and approximate results, |r_o(exact) — TD(approx)|/rD(exact), 
versus the absorption (i.e., detection) probability used to calculate To(approx) at the maximum. Same system as for the 
previous figure 

V. SOME EXTENSIONS 

In this section we present miscellaneous extensions of the previous formalism and results: for bound states, for 
dwelling into states rather than in a spatial region, and for multiparticle systems. 

A. The Harmonic oscillator and systems with bound states 

If we were to try the direct application of definition ([T]) to the case of the harmonic oscillator we would soon run 
into trouble due to divergent integrals. A natural option is to restrict the total time to one period of the system, as 
follows: 

Td = f T/2 dTe iST / n xn(x)e- iST / n , (41) 

J -T/2 

where the Hamiltonian is p 2 /2m + mw 2 i 2 /2, and T — 2ir/u>. We can thus write the dwell-time operator, restricted 
to one period, as 



oc 

Td =T^2(n\x D (x) \ n) \n) (n| 



(42) 



The eigenvalues, T (n \xd(x)\ n), can be understood as the period times the proportion of that period spent by the 
particle in the stationary state in the region D. This suggests an alternative useful quantity for systems whose 
Hamiltonian has a purely point spectrum. Define 

td = Hm i f T/2 dre^/ n XD (x) e" 1 ^ . (43) 

Let the Hamiltonian be given by 

H = Y J En\E n ,a){En,a\ , (44) 

with a a degeneracy index. Then it is easy to compute 

td = Yj \ E ^ a ) (E n ,a\xD{x) \E n ,P) (E n ,(3\ , (45) 

that is to say, the operator that gives us the fraction of time spent in region D is diagonal in the energy basis, up to 
degeneracy in energy. 



B. Times of Residence 



The construction above suggests an extension to the time of residence, which is an analogue of the dwell time valid 
for systems in which the concept of a region in space is not applicable. The question we now purport to answer is 
the following one: how much time has a system spent in a state \ip}7 Or, alternatively, from time to T, what is 
the proportion of time the system has spent in that state? Let then the projector associated with the state \if)) be 
denoted by P, and the unitary evolution operator by U(t). We define the time-of-residence operator as 

fv,(0;T) = / dtU\t)PU(t) . 



For example [47| , consider a two- level system with Hamiltonian 

q h ( Q 
H - 2 U 

and the state \tp) = ( „ J . The time of residence gets written as 



r, (0- T) ( 7 + m ^ (m g (1 - cos (OT)) 
T ^ U ' J j -^(1- cos (OT)) f-isin(nT) 

This entails that the measurement of this quantity would inevitably lead to one of the following two (eigen)values: 

T sin (CLT/2) 

T+ = — i — . 

± 2 n 

The fact that only two values can be obtained in each realization of the experiment, and not a continuous distribution 
of time of presence in the state = ( n ) , has been attributed to a predictive character of the measurement involved, 



v t 

as opposed to an observation that extends over time through the coupling with a weakly interacting clocking system 

As a matter of fact, the average time spent by a particle prepared at instant in a generic state ranges from 
t_ to t + for the eigenstate of P. Represent a generic pure or mixed state p by a vector in or on the Bloch 

sphere, p — (1 + r ■ c?)/2; then the average time of residence in state over ^ ne ti me interval [0,T] is 

T/2 + [(1 — cos (HT))r y + sin (flT) r x ] /2fi, which, as stated, ranges over [t_, r+]. 

It should be noticed that the operator of time of residence from instant to instant T is generically not stationary; 
for the example at hand we have 



H,t^(0;T) 



'&( 1 - cos (fiT) i sin (OT) 
2 I -isin(OT) -l + cos(nT) 



Notice that whenever T = 2nir/£l, with n a natural number, the corresponding time of residence does indeed commute 
with the Hamiltonian; this is in keeping with the expressions of Section [V Al 



C. Multi-particle systems 



The definition J!} admits a straightforward extension to systems of many particles, using the formalism of second 
quantization, in which we perform the substitution 

Xd(x) — > no = / dxa x a x , (46) 
Jd 

where a x is the annihilation operator at point x. For free particles, that is to say, for a Hamiltonian of the form 

H= d*— n*=/ d*— ajo*, (47) 



with at the operator that annihilates a particle with wavenumber fc, we can compute the dwell-time operator as 

(48) 

It is also feasible to write a dwell-time density operator, which for the free case reads 



i{a k + Ji eikL sin ( kL ) ai k a k 



6a(t) = / dk(^[8(t-t+(k)) + S(t-t-{k))]ala k 



oo 



[* (f - t+(*0) -«(*-*-(*))] 4<*-4 • 



It is easy to see that (i/)\IL a (t)\'ip) reduces to the dwell-time density for the one particle marginal wavefunction. The real 
multiparticlc aspect of this construction will only be accessible through dwell time - dwell-time correlation functions, 
for which a suitable operational interpretation needs to be built. 

In addition, we note that for multiparticle systems, the alternative question might be posed, as to what is the time 
during which a given number of particles n can be found in the region D. This naturally leads to the introduction of 
the dwell-time operator for n-particles 

/oo 
dtU\t)8n D ,nU(t) 
-oo 

= — dt d9UHt)e i9 ^ D ^U(t) , (49) 
2tt J-oo J q 

where is the density operator restricted to the region D, defined in Eq. (|4"6"| . Letting p be the density matrix 
describing the state of the system, it is convenient to introduce the characteristic function, 

F(0;t) = Tv[pUHt)e ieflD U(t)} , (50) 

whose Fourier transform is the atom number distribution in the region D (22j, 

P D (n, t) = ± e- in9 F(8; t)d0 . (51) 
27T Jo 

The meaning of Po(n,i) is precisely the probability for n particles to be found in the spatial domain D at time t. 

Knowledge of the atom-number distribution can be used to compute the average dwell time for different number of 
particles, namely, 

/>oo />oo 

(T D (n))= dtTY[p(t)Sr lD , n }= dtP D (n,t) . (52) 



VI. RELATION TO FLUX-FLUX CORRELATION FUNCTIONS 

This section follows [13] , and examines a link between the dwell time and flux-flux correlation functions (ffcf ) that 
have been considered mostly in chemical physics to define reaction rates for microcanonical or canonical ensembles 
[26[. The motivation for this exercise is to relate the dwell-time distribution, and not just the average value, to other 
observables. 



A. Stationary flux-flux correlation function 



Pollak and Miller [40( have shown a connection between the average stationary dwell time and the first moment of 
a fief. They define a quantum microcanonical ffcf Cpm{t, k) = Tr{Re Cpm (t, k)} by means of the operator 

C PM (r,k) = 2irh[J(x 2 ,T)J{x 1 ,0) + J( Xi ,t)J(x 2 ,0) 

- J(x 1 ,t)J{x 1 ,0)-J(x 2 ,t)J(x 2 ,0)}S(E-H), (53) 

where J(x, t) = e lHt / n ^{pS(x — x) + S(x — x)p\e~~ lHt / n is the quantum mechanical flux operator in the Heisenberg 
picture, and p and x are the momentum and position operators. 

This definition is motivated from classical mechanics: Eq. (|53p counts flux correlations of particles entering D 
through xi (x 2 ) and leaving it through x 2 (x±) a time r later. Moreover, particles may be reflected and may leave 
the region D through its entrance point. This is described by the last two terms, where the minus sign compensates 
for the change of sign of a back-moving flux. Note that these negative terms lead to a self-correlation contribution 
that diverges for t — > 0. 

We shall first derive the average correlation time and show its equivalence with the average dwell time. We shall 
only consider positive incident momenta, and define Cp M by substituting 5(E — H) by d + (E — H) := 8{E — H)A + , 
where A + is the projector onto the subspace of eigenstates of H with positive momentum incidence. By means of the 
continuity equation, 

-— J(x,t) = —p(x,t) , (54) 
where p(x, t) — e lHt f h 8(x — x)e~' Ht ^ h is the (Heisenberg) density operator, Cp M (r, k) can be written as 

C^ M (r,k) = -2nh (^Xd(t)) _ 5 + (E - H) . (55) 



By a partial integration and using the Heisenberg equation of motion the first moment of the Pollak-Miller correlation 
function is given by 

Tr{^°°dTTC+ M (T,fc)J =Tr{2irhJ™dTXD(T)±[x D (0),H]8+(E-H)\ . (56) 

Boundary terms of the form limT-^oo t 7 £c(t), 7 = 0, 1,2, are omitted here and in the following. The contribution 
of these terms should vanish when an integration over stationary wavefunctions is performed to account for the 
wavepacket dynamics, as it is done explicitly in the next section. For potential scattering the probability density 
decays generically as r~ 3 , which assures a finite dwell-time average, but for free motion it decays as r _1 [3J|, making 
td infinite, unless the momentum wave function vanishes at k = sufficiently fast as k tends to zero (Tl| . as we have 
discussed before. 

Writing the commutator explicitly and using the cyclic property of the trace gives 

Tr{^°°dTTg+ M (T,fc)} 

= Tr 1 2wh J™ dr(--^xn(T)\xD (0)S+(E-H)\ , (57) 

and integration over r yields the final result, 

Trj^ drTC+ M (T,fc)| =2irKTr{x D (0)S+(E-H)} =J kk . (58) 

Expressing the trace in the basis \4>k) gives back the stationary dwell time of Eq. ©, i.e. the diagonal element of the 
on-the-energy-shell dwell-time operator, Tkk- 

The calculation of the average in |40l | is different in some respects: The coordinates x\ and x 2 are taken to minus 
and plus infinity, but it can be carried out for finite values modifying Eq. (8) of [iol ] accordingly; Formally there are 
no explicit boundary terms at infinity but a regularization is required in Eq. (16) of [40| . which is justified for wave 
packets; (c) 5(E — H) is used instead of 8 + (E — H). That simply provides an additional contribution for negative 



momenta parallel to the one obtained here for positive momenta; (d) In our derivation the average correlation time 
is found to be real directly, even though C~k M (j, k) is not selfadjoint, whereas in the real part is taken. (The 
discussion of the imaginary time average in [401 ] is based on a modified version of Eq. ([55]) .) 

Next, we will show that the second moment of the Pollak- Miller ffcf equals the second moment of T. This was not 
observed in Ref. (40| . Proceeding in a similar way as above, we start with 



Z = Ti- 



ll 



d^C+ M (T,£0 



(59) 



Integrating by parts twice, neglecting the term at infinity, using Heisenberg's equation of motion, and the fact that 
<pk is an eigenstate of H, the real part is 



J + J* 27TTO 



Tik 



dr(<fe | [xd (t)xd (0) + xd (0)*d (t)] \<p k ) 



(60) 



Introducing resolutions of the identity, 



ReZ 



hk J Q 



OO pOQ 



— oo J X\ 



dr\ dk' dx dx'e^-^^cPU^k'ixWvix'^kix'y. 



+ c.c. 



(61) 



where c.c means complex conjugate. Making the changes r — * — r and x,x' — > x',x in the c.c-term, it takes the same 
form as the first one, but with the time integral from — oo to 0. Adding the two terms, the r-integral provides an 
energy delta function that can be separated into two deltas which select k' = ±fe to arrive at 



Tr <^ Rc 



dTT 2 C+ M (T,k)\ = 



<in 2 m 2 



h z k 2 



x 2 s 2 

dx |^fe(x)| 2 



dx <j>l(x)</>-k(x) 



(62) 



In other words, the relation between dwell times and flux-flux correlation functions goes beyond average values and 
Cp M (r, k) includes quantum features of the dwell time: note that the first summand in Eq. (|62p is nothing but 
(Tfcfc) 2 , whereas the second summand is positive, which allows for a non-zero on-the-energy shell dwell-time variance 
(T 2 )fcfc — (T^) 2 . We insist that the stationary state considered has positive momentum, 4»k(x), k > 0, but this second 
term implies the degenerate partner <fi-k{x) as well, and is generically non-zero. 

We shall see in the next section with a more general approach, that these connections do not hold for higher 
moments. 



B. Time-dependent flux-flux correlation function 



A time-dependent version of the above flux-flux correlation function can be defined in terms of the operator 

/oo 
dt [j(x 2 ,t + T)J(x 1 ,t) + J(x 1 ,t + T)J(x 2 ,t) 
-oo 

- J(xi,t + r)J(x 1 ,t) - J(x 2 ,t + T)J(x 2 ,t)] , (63) 



which leads to the flux-flux correlation function 



C(r) = (ReC(r)) 



(64) 



where the real part is taken to symmetrize the non-self adjoint operator C(r) as before. 

As in the stationary case, Eq. (|64p counts flux correlations of particles entering D through x\ or x 2 at a time t and 
leaving it either through x\ or x 2 a time r later. Moreover we integrate over the entrance time t. It can be shown 
that the first moment of the classical version of Eq. (|63| . where J is replaced by the classical dynamical variable of 
the flux, gives the average of the classical dwell time. 



As in Eq. (|55p . we may rewrite C(r) as 



f°° d d 
C(t)=-J &t—XD{x,t + T)—x D {x,i). (65) 

From here we note that the ffcf C(r) is not normalized, 

f°° d 

dr C(r) = J dt xd (t) —x D (t) = 0, (66) 

as a result of the self-correlation. 

Next we derive the average of the time-dependent correlation function. With a partial integration one finds 

/oo poo ^ 

dt dTx D (x,t + T)—x D (x,t). 
-oo JO at 

A second partial integration with respect to t, replacing d/dt by d/dr and integrating over r gives 

/oo />oo 
df xl @, t)= dt xd (x, t)=f D , (67) 
-oo J —oo 

where Xjj — Xd has been used. Eq. ([6"T|) generalizes the result of Pollak and Miller to time-dependent dwell times. 

A similar calculation can be performed for the second moment of C(r). After three partial integrations with 
vanishing boundary contributions to get rid of the factor r 2 one obtains 



y* OO pOO />oo 

/ dTT 2 d(r) = 2 dt dTX D (x,t + T)xn(x,t) . 

JO J -oo JO 



10 J-oo JO 

Making the substitutions t + r — > t and r — » — t in the complex conjugated term, we find 



Re / dTT 2 C(T) =T 2 . (68) 



o 

C. Example: free motion 

For a stationary flux of freely-moving particles with energy Ek, k > 0, described by 4>k{x) = (x\k) = (27r)~ 1 / 2 e l/cx , 
the first three moments of the ideal dwell-time distribution on the energy shell are given by 

_ mL , . 

T fefe = , (69) 

2 _ m 2 L 2 f sin 2 (fcL) \ 
/t3\ _ m 3 L 3 f sin 2 (fcL) \ 

As proved above, the first two moments agree with the corresponding moments of the Pollak-Miller ffcf, but for the 
third moment we obtain instead 



Trine rdTT 3 C+ irk)} - fl 3[1 ± ^ kL)] I 3sin(2fc£) 



(72) 



In Fig. [K] the first three moments are compared. The agreement between (T 3 )kk and Eq. |72p is very good for large 
values of fc, but they clearly differ for small k. Nevertheless, the agreement of the first two moments suggests a similar 
behavior of II(r) and C(r). 

To calculate H^(t) for a wavefunction ip(k) '■= (k\ip) with only positive momentum components we use Eq. (jlip 
and the explicit forms of the eigenvectors, Eq. (1141) and eigenvalues t±, Eq. (fT5|) . 

n " (r) "2^ 7 ? ± ^OT' (73) 
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FIG. 5: Comparison of the first three moments: 



(T 2 )kk and (T 3 )^ (dotted-dashed line) with the corresponding moments of 



the fiux-fiux correlation function, for a free-motion stationary state with fixed k. (Tkk) 2 is also shown (dotted line), h 
and L — 2 



where the j-sum is over the solutions fej(r) of the equation Fj(k) = i 7 (fc) — r = and the derivative is with respect 
to k. 

We use the following wavefunction [29(, 

${k) = Af(l - e -"k 2 } e -(k-k y/l4{Ak) 2 ] e -ik XoQ ^ ; (74) 

where Af is the normalization constant and Q(k) the step function. For the free flux- flux correlation function we write 



/>oo />oo 

C(r)=Re/ dk dk' ip*(k)ip(k')(k\d(T)\k') 
Jo Jo 



and CW(t) = (k\C(T)\k'} in the free case 



where 



jn d^ 

— —6(k — k')—— 2g(hkr lm) — gthkr/m — L) — gikkr/m + L) 
2-Khk <1t 



1/2 



(75) 



(76) 



(77) 



The result is shown in Fig. [6] The ffcf shows a hump around the mean dwell time but it oscillates for small r and 
diverges for r — > 0. As discussed above, this is due to the self-correlation contribution of wavepackets which are at 
x\ or xi at the times t and t + r without changing the direction of motion in between. A similar feature has been 
observed in a traversal-time distribution derived by means of a path integral approach [lil ]. 

In contrast, n^(r) behaves regularly for r — > 0, but shows peaks in the region of the hump. This is because the 
denominator of Eq. ([75)) becomes zero if the slope of the eigenvalues t±(k) is zero, which occurs at every crossing of 
t+{k) and t-{k). 

The distribution 7t^(t), see Eq. (|?7|). is also computed: it agrees with C(r) in the region near the average dwell 
time and it tends to zero for r — > 0. However, it does not show the resonance peaks of 11^, (r). 

In absence of a direct dwell-time measurement, the physical significance Tjj and to depends on their relation to 
other observables. The present results indicate that the second moment of the fiux-fiux correlation function is related 
to the former and not to the later, providing indirect support for the physical relevance of the dwell-time resonance 
peaks, but other observables could behave differently 



D. Approximations 

While the previous results bring dwell-time information closer to experimental realization, the difficulty is translated 
onto the measurement of the ffcf, not necessarily an easy task. A simple approximation is to substitute the expectation 
of the product of two flux operators by the product of their expectation values, which are the current densities. Using 




FIG. 6: Comparison of dwell-time distribution II(r) (dashed line) and flux- flux correlation function C(r) (solid line) for the 
freely moving wave packet (|74p . Furthermore, the alternative free-motion dwell-time distribution 7t(t), Eq. (127[l . is plotted 
(circles). We set % = m = 1 and |a;o| large enough to avoid overlap of the initial state with the space region D = [0, 45]. fco = 2, 
Ak = 0.4, and a = 0.5 




FIG. 7: Comparison of the relative error of (To) using the approximation Co(t) instead of C(t) for free motion, a = 0.5, 
h = m = l and L = 100 



the wave packet of Eq. (fT4"|) . we have compared the times obtained with the full expression (frj4"|) and with this 
approximation in Fig. [7] The two results approach each other as — > 0, also by increasing L and/or /co- 

The exact result can be approached systematically, still making use of ordinary current densities, as follows: First 
decompose J(xi,t + r)J(xj,t) by means of 



1 = P + Q 

P = |V>>M : 



so that 



J(x h t + T)J(xj,t) = J{xi,t + T){\il)){rj}\ +Q)J(xj,t) . 
It is useful to decompose Q further in terms of a basis of states orthogonal to and to each other, 



(78) 
(79) 

(80) 
(81) 



that could be generated by means of a Gram-Schmidt orthogonalization procedure. Now we can split Eq. ()63|) . 

C(t) = Co(t) + C 1 (t) , (82) 

where Cq(t) has the structure of C, but with P inserted between the two flux operators in each of the four terms. 
Similarly G\ (r) has Q inserted and can be itself decomposed using Eq. (|8l"|) . 

We define C(r) = Cq(t) + C\{t) by taking the real real part of (i/j\Cq(t) + Ci(t)\i1j). Co is the zeroth order 
approximation discussed before and only involves ordinary, measurable current densities [9( . The non-diagonal terms 



from C\, (ip\J(xi,t)\tjj^)(tp l ^\J(xj,t + r)\if>) can also be related to diagonal elements of J by means of the auxiliary 



states 



since 



IV-i) = \i>) + \^f) > 
\*h) = |V) + #?>, 

1^3) = H)-^f), (83) 



(^|J(x,t)|^ Q ) = ±{ik.\J(x,t)\1> l ) - ±(ik\J(x,t)\ik) 

- ±{fo\J(x,t)\fo) + fa\J(x,t)\ifa) 

- l(to\J{x,t)\<k) . (84) 



VII. FINAL COMMENTS 



The quantum dwell-time distribution of a particle in a spatial region, and its second moment present non-classical 
features even for the simplest case of a freely moving particle, such as bimodality (due to two different eigenvalues 
for the same energy) with the strongest deviations from classical behavior occurring for de Broglie wavelengths of the 
order or larger than the region width. Progress in ultracold atom manipulation makes plausible the observation of 
these effects, and motivates further effort to achieve an elusive direct measurement of the dwell times, or to link the 
distribution and its moments to other observables. We have in this regard pointed out that the flux-flux correlation 
function provides access to the second moment. The potential impact in cold atom time- frequency metrology [Il,[44[, 
and other fields in which the dwell time plays a prominent role (such as conductivity pj], or chaos) remains an open 
question. 
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